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Abstract 

We extend the two particle theory of disordered systems within the coherent 
potential approximation CPA to obtain weighted contributions to averaged 
two particle resolvents which arise from separate alloy components. Starting 
from first principles in a model of diagonal disorder and the single site approx- 
imation for a binary substitutional alloy A C B\- C we extend the approach of a 
fundamental paper by Velicky to evaluate various weighted forms of a general 
class of two particle Green's functions. Applications in a wide range of linear 
response theory are discussed in detail as well as the behavior of the weighted 
functions in a strong disorder limit. To exemplify our analytic calculations 
the optical absorption in a disordered model-alloy is studied numerically. 

PACS numbers: 71.10,+x 78.50.-w 72.10.Fk 71.55.-i 



1 



Typeset using REVTgX 



I. INTRODUCTION 



In recent years the understanding of the effects of disorder in the physics of metals, 
semiconductors and amorphous systems has made a tremendous progress. This vigorous 
development was motivated to a great extent by a thorough understanding of how strongly 
disorder effects determine the behavior of real physical systems. The success of the various 
analytical descriptions which have been considered, however, has always been determined 
by the relative simplicity to which approximations could be reduced, in order to keep the 
theory tractable analytically, without losing its capability to account for the most important 
physical aspects. One of the most successful approximations to match these requirements 
has been the coherent potential approximation CPA developed by Soven0, TayloiH, and 
extended by LeatlM Velickyi~i and many others (see Ref. for a review). 

In contrast to many other approximations for the procedure of configurational averaging 
in disordered systems, the CPA is capable of interpolating correctly between the limits of 
weak0 and strong disorder!! as well as low and high impurity concentrations. Therefore 
it is also able to predict accurately the formation of impurity bound states turning into 
split-off impurity bands as the impurity concentration is increased. Extensions of the single 
site CPA to include scattering from clusters of impurities are still today the only existing 
analytic theories that allow for a calculation of the density of states in disordered systems 
such that a reasonable versimilitude is attained^. The main difficulty with the CPA, is the 
relative complexity of the self consistent equations which have to be solved in more accurate 
extensions of the theory. 

In the present paper we propose a more differentiated analysis of two particle properties 
within the CPA building on the single site approximation (SSA) in which it was developed 
originally. Our novel treatment should be especially useful for applications to binary substi- 
tutional alloys in all regimes of disorder. Following an earlier paper by Velickyi, in which the 
regular two particle theory within the CPA was first introduced and which will be referred 
to hereafter as I, the theory of weighted single particle resolvents^ is extended to a properly 
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weighted two particle theory. Weighting of a Green's function in this context means that 
through the application of appropriate operators to the unaveraged resolvent, restrictions 
are made on the type of alloy component of either or both of the sites on which the particle 
starts and terminates its motion. Upon averaging, this results in a statistical weight being 
attributed to the averaged unrestricted resolvent with possibly separate additive terms. 

On the average, therefore, contributions to a full two particle Green's function in terms 
of constituent components can be resolved. This allows for a better understanding of how 
several physical processes contribute to a cumulative behavior as the principal parameters 
of the system are varied and it can therefore be used in an analysis of further effects which 
differentiates between these components. 

We have structured the paper as follows: In section II the important features of the single 
particle CPA for diagonal disorder are recapitulated and relations which are important to the 
calculation of the corresponding two particle Green's functions are established. In section 

III we calculate the Fourier transforms of a class of weighted two particle functions which 
are kept as general as possible. For that reason only a representative choice of weightings are 
calculated explicitly, since other weights can be obtained in an analogous fashion. Section 

IV is devoted to possible applications of weighted two particle functions in linear response 
theory, and the peculiarities of two different classes of such functions are discussed in detail, 
which we have selected to cover a large range of conceivable applications. 

In the third part of section IV the behavior of one class of functions discussed before is 
examined in a split band limit (strong disorder). Section V is devoted to a numerical study 
of the splitting into several components of an interband absorption spectrum of a disordered 
alloy. Section VI in conclusion discusses the implications and possible further applications 
of the results obtained throughout this work. 
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II. SINGLE PARTICLE PROPERTIES 



We consider a binary substitutional alloy A c Bi_ c on a simple totally disordered 
monoatomic lattice on which each site is occupied by either an A— or a .B-component with 
probabilities c and 1 — c, respectively. Since by symmetry it is only necessary to consider 
concentrations between < c < 0.5, we restrict our investigation to this region and define 
A as the impurity and B as the host component of the alloy. 

To simplify further calculations we consider diagonal disorder only, i.e. a Hamiltonian 
of the form 

H = H + U = H + J2U n (1) 

n 

where H represents the periodic part and the U n are single site contributions to a random 
disorder potential which assume either the value = V or f/jf = with probabilities c 
and 1 — c, depending on whether the site n is an impurity or part of the host material. 
In principle one could also have introduced a symmetric model for the disorder, but as it 
turns out the amount of algebra is somewhat reduced by the asymmetric definition, while 
switching from one form to the other does not present any difficulty 

The propagator of the disordered medium G relates to the one of a pure host medium 
(pure B-phase) g through the Dyson equation 

G = g + gUG = g + gUg + gUGUg. (2) 

The CPA for a disordered medium is introduced by the usual method^H of placing the 
impurities in a self-consistent medium such that the the averaged propagator of this effective 
medium fulfills the relation 

G = g + 9 m (3) 

where S is the CPA single particle self energy. S itself is determined by the self consistency 
condition that the average of the total scattering of a particle in the effective medium be 
zero. This total scattering is described by the equation 
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G = G + GTG. (4) 

where T is corresponding T-matrix of the problem and the self consistency condition is 
therefore (T) = 0. In the SSA an additional requirement is made through the condition 
that the total scattering off a single site n be zero. This scattering is described by the single 
site contribution to the T-matrix T n which is defined as 

T n = (U n - £„) [1 + FT n ] (5) 

where F is the site diagonal average propagator and U n and S n are the single site decom- 
positions of U and E, such that U = J2 n U n and S n = J2 n ^n- The average of the disorder 
potential U alone amounts to (U) = cV This along with the average of (|5]) set to zero 
determines E to be 

(6) 



1 - (V - £)F' 

One can now define conditional or weighted propagators^ which explicitly describe the 
propagation of a particle between partly or completely specified types of sites by multiplying 
a normalized version of the random potential onto them. For example, 

G " = "~v (7) 

describes the motion of a particle commencing on an impurity site and ending at an arbitrary 
other site in the medium, since U will be zero if the first site of the function is a host. 
Similarly, 

G" = ^ (8) 

describes a situation where both sites are required to be impurities for the function not to be 
zero. Upon averaging over all configurations, the Green's functions become translationally 
invariant and the following relationships between the averaged weighted and unweighted 
functions can be obtained 
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V 



G' 



E-cV 



(10) 



The second term corrects the site diagonal elements when G n = G % and uses identity @. 
From here on, other weighted functions can be calculated by probability conservation. It is 
found that 

E' 



G h 

Qih 
£<hh 



G-G l = 



G 

E / E , 
1 — — ) G 



V 



V 

E iV<5 



E- C y 



(ii) 

(12) 
(13) 



It is our goal in the present paper to establish the two particle analogues of these weighted 
Green's functions, i.e. jointly averaged products of such functions including the coherent 
scattering which induces correlations in the joint propagation of two particles. 

In order to deduce these jointly averaged weighted functions, we first need to obtain the 
non- averaged weighted functions, starting with G l and G u , in a representation such that no 
products between the disorder potential U and the unweighted single particle function G 
occur. To accomplish this, one can simply employ equation (fj) which yields 

g-'G-i 



G 



V 



riii 1 



and by means of (f|) these go over to 

1 



g-'Gg- 1 - g 



-i 



U 



(14) 
(15) 
(16) 



1 



g~ l G + g~ x GTG — 1 



9~ X Gg~ x + g~ l GTGg~ L - g 



-i „-i 



U 



(17) 
(18) 



A further single particle identity resulting from (||) which will prove to be very useful is is 



(19) 



6 



III. TWO PARTICLE THEORY 



In this section a general weighted two particle theory involving two different bands is 
established from which the case of two particles moving in a single band can also be imme- 
diately obtained. The calculation for the unweighted two particle function has been made 
in I. As a main result of that work the appropriate vertex corrections for the CPA were 
obtained which account for the coherent scattering processes of two particles that arise in 
the otherwise non-interacting two particle function through the averaging process. 

We follow the outline of Velicky's reasoning to obtain the proper weights for two choices 
of the functions (G^CG^ v ). The labels /i, v, //, v' G {i, h, 0} indicate the kind of weight - 
either an impurity a host or no weight - which is attributed to the first and second site of 
the respective function. C represents a generalized operator coupling the two single particle 
functions and a and b label two possibly different bands on which the respective single particle 
resolvents are defined. The positions of the weights within one single Green's function are 
thereby important, since the disordered medium before averaging is neither homogeneous 
nor isotropic and thus the non-averaged Green's functions depend non-trivially on both 
arguments. In principle there now would be 80 (!) possible ways of applying specific weights 
to these functions before averaging. However, in most cases, even if different bands are 
involved, only two particle functions having an equal kind of weighting on its single particle 
constituents will be needed in most applications. Later, we will consider two particular 
examples of two particle functions which find frequent use in linear response theory where 
the operator C is diagonal in real space and also only diagonal elements (or sums of diagonal 
elements) of the functions defined above are used. For both cases considered, the total choices 
of weightings reduce to only five different ones, since the first and the second site of the first 
function will be the same as the second and first site of the second function, respectively, 
which implies that they must also pairwise bear the same weighting label. 

As an example we calculate for a most general choice of C only the two functions 
(G™CGfr) and (GfCG® 1 ) in this section, since they will prove to be the most useful types of 
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weightings for the cases discussed thereafter and all other ones could be obtained in complete 
analogy to the calculations presented here. 

Using the identities ( |TTD and ([18]), the single weighted two particle function can be written 
in terms of the single particle functions as 

(GfCGf) = V-X'Wg-'Ga + g-'GaTaGa -1]C \G b T b G b g b l + G b g^ - ll ) (20) 



and the double weighted one as 

(G^CGt) = V^V b -*([{g-'G a g- x - g' 1 ) + g^G a T a G a g^ - U a 
x \g b - l G b g b - 1 - g b l ) + g^G b T b G b g^ - U b \). 



C 



(21) 



We calculate the double weighted function first, since it provides the more difficult task and 
from its solution it is straightforward to derive the one for the single weighted function as 



well. The problem of evaluating ([H]) is divided into two parts, the first one involving all 
terms not containing the matrix U, and the second one containing the remainder, i.e. 



V a V b \GaCG") = JC + M 



(22) 



K 



g-'Gag- 1 - g~ x ] C Ig^G.g, 1 - g b l ] + g^Ga^G^Gg^G^G^ 1 (23) 



M = (U a CU b ) + [g- l C(JJ b ) - g-'Gag-'CiU,) - g^G a (T a G a g^CU b ) + (a <- b)} (24) 

where in /C the terms involving an average over a single T-matrix have vanished, which is 
the standard CPA condition and (a <-> b) indicates that the labels are exchanged and the 
corresponding expressions reflected around the operator C. We evaluate /C first since it is 
the term needed in the wider range of applications. 
With the identity (|19D we find that 

g^G a g- 1 - g- 1 ] C [g b l G b g b l - g b l ] = E^g^GaCG^ 1 . (25) 

The T-matrix can be decomposed into its single site contributions T n as 



T = J2T n +J2 T nGT m + E T n GT m GTi + .. 



(26) 



Thereby the characteristic exclusions in the sums prevent the particle from scattering twice 
in sequence on the same site and T n satisfies equation (||). 

As shown by Velicky, T can then be replaced in two ways by a closed set of equations, 
namely 

T = J2Qn = J2Qn (27) 

n n 

where 

Q n = T n (l + Gj2Qm) (28) 

and 

Qn = (1 + E QmG)T n . (29) 

Due to the requirement that (T) = and the single site decomposition of T from (|26| ) it is 
possible to decompose averages on different sites to give 

= <r> = = £(T n > ( 1 + G £ (Q m > ) (30) 

This also implies that (Q n ) = (Q n ) = 0. A vertex function T can now be defined similar to 
I such that 

K, = g- l G a {Y,j: b C + T)G b g? (31) 

where now 

r=(T a G a g~ 1 Cg b - 1 G b T b ). (32) 



T can be manipulated along the lines of I by using (|27D to yield 

r = Y,T,(QnG a 9a 1 C 9b - 1 G b Q b m ). (33) 

n m 

By means of ( p8|) - (^) this can be cast into the form 

T n = (T«G a g-\C + g a J2 T p g b )g b l G h T h n ) (34) 
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where now T = ^2 n T n , since from the decoupling introduced in (|30|) one gets 
(Q a n G a g~ l Cg^G b Q b m ) = (Q a n G a g- l Cg^G b Q b n )5 n>m . The only difference to the correspond- 
ing expression in I, (cf. (47) there) is that J2 P ^ n T P is surrounded by the propagators of the 
pure medium g a / b here. We can then use ( pT|) and (|34| ) to obtain 

Tn = (TnG a g a 1 Cg b 1 G b T b ) + (T^g a K,g b T b ) — Y la Y lb (T%G a CG b T b ) — (T%G a T n G b T b ). (35) 

At this stage we are able to find T n and therefore also K,. As we cast our model into a 
site representation we obtain T n —\ n)t n (n | and hence 



T n =| n)A{n | [g a fcgb + G a g a Cg b G b - Ej: b G a CG b \ | n)(n 

where A is the irreducible vertex part derived by Velicky 

(t«( Zl )t b n (z 2 )} 



(36) 



(37) 



The vertex A can be regarded as being intrinsic to the CPA, since it does not depend on 
the particular form of the operator C. Substituting (^) into (|3T| ) yields 

= S a £ b g a 1 G a CG b g b 1 + 
+ kg- l G a ]T | n)(n \ [g a )Cg b + G^Cg^Gy - Z a Z b G a CG b ] \n)(n\ G b g b \ (38) 

n 

Multiplying by g a and g b from the left and right, we solve for the diagonal elements (m | 
gjigb | m) 

(m | g a fcgb | rn) = S a E 6 (m | G a CG b \ m)+ 



+ A ^2(m | G a | n) (n \ g a JCg b + G a g a 1 Cg b 1 G b - T, a T, b G a CG b \n)(n\G b \ m) (39) 

and hence also solve fl38|) . At this point it is helpful to visualize the form of the operator C 
in the site representation, which in the most general case can be written as 



l,m 
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(40) 



For convenience we introduce the short notation 



T 

•J n. 



(n | G | m), 
(n | g" G | m) 



(41) 
(42) 



where Fq = F as already defined in (||). It should be noted that F n _ m and T n -m are of 
different dimensions and their definition has been chosen to reduce the algebra as much as 
possible. Furthermore, these expressions show that equation fl39[) only contains translation- 
ally invariant quantities and hence it can be solved by Fourier transformation. The following 
Fourier transforms are introduced 

(43) 

(44) 
(45) 
(46) 
(47) 
(48) 

The units of ak,ak,bk,Ak are unity, that of Ak is J -2 and that of Ck is J 2 . Inserting into 
© yields 

a k + AA k (a k - a k ) 





m 


ikRm {m | G a CG b | 


m) 


Oik 


m 


m | g^GaCGbg^ 1 


m) = 




^ ^ g ikRm 


S a S b (m G a CG b 






m 




n 


A k 


_ -ikRm fpa -pb 
/ j c m — m 
m 




A k 


_ ST^ —ikRm <pa >pb 
/ j e • r m- r —m 

m 




h 


_ ^ e -ikR m ^ 
m 


m | | m) 




Ck 


_ e ~ikR. m ^ 
m 


m | K | m). 





l - A A- 

and hence by means of ( |38| ) 

S a S 6 + A(A k - EaE b A k ) 



(49) 



Cfc = oi k - 



1-AAj 



S a S fe + AD k 



1 - A A, 



(50) 



where 



F>k — Ak — Xa£b^4* 



e -ikR n 

III 



1mm ~\~ S a r )nmF m _ n + Sfe 1mnF n _ 
n n 



(51) 
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Similarly, one can evaluate the term M. from (p4[), but it turns out that little simplification 
can be made until an explicit form of C is known and therefore we postpone its evaluation 
to the next section. 

For now, we proceed to calculate the single weighted function (G^CG® 1 ). From ( f20|) we 

get 

VaVbiG^CG® 1 ) = g a 1 G a CGi > g b 1 + C — g a 1 G a C — CGbg b 1 + g a 1 G a {T a G a C G\ } T h )Gi ) g h l . 

(52) 

Using (|T9|), this can be recast into 

V a V b (G%CG?) = E a E b G a CG b + g^G ^0^0^)0^ (53) 

which can be readily solved since the term (T a G a CGbTb) exactly corresponds to the vertex 
part of the unweighted function which is known from I 

(T a G a CG b T b ) = J2 I n)A(n | (G a CG b ) \n)(n\. (54) 

n 

Thus, the Fourier transform 

4 = V a V b £ e ikR - (m | (GfCG™) \ m) (55) 

m 

can be obtained using the solution for the unweighted two particle function and equations 
(H), (H) and (H) to yield 

4 = a fc + AA SaSb( ;i AAfc) . (56) 
Using fl5T] ) this can be written as 

afc T> a T> b + AD k 
4 " l-AA k ■ (5?) 

It is important to indicate at this point that the solution for a single weighted function in 
which the weights have been swapped to the inside, i.e. (G^CGl®), will not be the same as 
the one just obtained. A similar type of calculation yields instead of ( j56f) 
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d' k = a k + AA k l ® k AAk , (58) 

where d' k = Y, n e ikRm (m \ (G^CGf) \ m). 

Since we have now obtained expressions for the Fourier transformed site diagonal el- 
ements of some of the impurity weighted two particle functions, we are also able to find 
their respective off diagonal parts by inserting into the appropriate Bethe-Salpeter equa- 
tions ((|3"5D in the case of the double weighted function). Other types of weighted functions 
can be calculated from here by setting up the corresponding equations in analogy to ( 20 ) 
and ( PH ) from the single particle functions. In particular, to obtain the corresponding host 
weighted functions one can apply the weights (1 — U a /b/V a /b) to the single particle resolvents 
from either side, but the form of this weight already suggests that host weighted properties 
can always be deduced by adding and/or subtracting corresponding unweighted/impurity 
weighted functions. For the most general case the set of closed equations will be quite large 
and for that reason we will in the following deduce these equations only in some more specific 
cases. 



IV. APPLICATIONS IN LINEAR RESPONSE THEORY 

The most obvious application of a two particle theory as introduced above is in linear 
response theory, i.e. in calculating weighted susceptibilities. The study of weighted functions 
in this context helps to determine how such quantities as susceptibilities and transport 
coefficients are constituted (on the average) from processes on different components of the 
alloy. A further advantage of such a differentiation is that it allows for a more refined 
treatment of further renormalizations to the considered quantities once further interactions 
are introduced into the problem. The standard expression employed in linear response theory 
is a generalized Kubo formula 

XcM,cw(*) = / ^d V S Q (z^, V )Z- 1 Tr b {e^ H C^(5^^ - H a )C^5^( V - H b ))}. (59) 
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Thereby XcW c< 2 ) stands for a generalized type of susceptibility characterizing the linear 
response of an observable to an external perturbation coupling into the Hamiltonian 
through the operator C^. So denotes the zero order associated two particle Matsubara 
function in the Lehmann representation^ and 5^(£ — H a /b) is the corresponding spectral 
function of the full resolvent G^j h {^) = H^,J£ — Ha/f,] -1 ^^, defined on the subspace (band) 
a/b, with the weighting operators E^, b applied to it which assume the values S^ 6 = U a /b/V a /b, 
^a/b = (1 — U a /b/V a /b) and H°, 6 = 1. Tr b denotes the trace over the subspace (band) b. Our 
particular interest focuses on the operators C' 1 ' and which we now take to have the 
following diagonal form 

C (1)/(2) = E7i 1)/{2) (^^)l^)(^l- (60) 

m 

To examine the behavior of the weighted functions further under such a constraint it is 
useful to consider the two cases: 

7 m = const x 5 m>r (61) 

where r is an explicitly specified site and 

7 m = const. (62) 

The first case arises in a treatment of the response from local interactions such as for example 
the space dependent exchange interaction between spin impurities (RKKY interactional) 
submerged in a system otherwise containing non-magnetic disorder^. The second case is 
widely used in calculations of the optical response in metals and semiconductorsSEiil, since 
the exciting optical fields can be taken as uniform over space with approximate momentum 
q = and the optical matrix elements describing transitions between bands of different 
angular momentum symmetry are usually well approximated as constants. 

In both cases the restriction to only consider diagonal elements of two particle functions 
reduces the set of closed equations for all possible weightings to 

(G h a CG h b ) = (G a CG b ) - (GlCGl) (63) 
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(Gf CG\ h ) = <G*CG*> = (GlCGl) - (G™CGf) 



(64) 



/ fihh ft /^ihh\ / r~ih fi r~ih\ / fiihfifiih\ 



(65) 



It follows, that for this case the calculation of only (G"CGf) and (G l a CG % b ) is sufficient to 
also obtain the remaining mixed and host weighted functions. We start to consider the case 
of (|T]) first. 



The response of a system at some point r\ in space to a perturbation applied at point 
r 2 has a wide range of applicability. Disregarding the coupling constants which simply scale 
the result we assume that =| 2) (2 | and =| |, where the choice of the origin 
is arbitrary. In a slightly different notation the calculation of the double weighted function 
corresponds to evaluating (G"(ri, r 2 ; t)G^(r 2 , T\\ £')). 

For the RKKY interaction mentioned before, which couples spins at a distance T\ — r 2 
through the electronic spin susceptibility of electrons in the conduction band, the two- 
particle time dependent function is to be Fourier transformed and to be taken at two iden- 
tical single particle energieslll. If the spins are at impurity sites, the weighted functions 
must be used to describe the problem adequately. Earlier treatments^ neglected the ver- 
tex corrections and it is proposed to investigate the effect of their inclusion in a separate 
publication. 

In this as in other problems the self-interaction r% = r 2 can be neglected for most 
purposes, which implies that it will be sufficient to calculate only the term K, from ([23|), 
since all terms occurring in Ai from ( p4|) will vanish for r% ^ r 2 , since the disorder potential 
U is a diagonal matrix which vanishes identically for off-diagonal terms. To clarify the 
meaning of the quantity which is obtained through this special choice of and C^ 2 \ we 
look at how the corresponding Bethe-Salpeter (B-S) equation can be rewritten for the case 



A. Electronic Susceptibility 
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of a unweighted two particle Green's function as introduced in equations (22) of Ref. [4[ In 
the single site approximation the B-S equation for the unweighted function reads 

(G„(l, 2)G 6 (1, 2)) = (G$(l, 1; 2, 2)) = (G„(l, 2))(G 6 (1, 2)) + 

+ A £(G a (l, n )> (G 6 (l, n))(Gi? (n, n; 2, 2)) (66) 

where we have rewritten 

(2 | G | | G b | 2) ee G 8 (1,2)G 6 (1,2) (67) 

and 

(2|(G cl |l)(l|G b )|2) = (Gi 2 b ) (l,l;2,2)). (68) 



For the Fourier transforms from before (E3h - (46) this amounts to 



a k = S a S 6 e- lkR ™G a (m)G b (m) (69) 

m 

A (70) 



m 



m 

= E e ~ tkRm [X a X b G a (m)G b (m) + 5^(1 + S a F a + E b F b )] = A k . (71) 

m 

Using flTU|) and (|7T| ) in connection with ( pB[ ) and floT^ ) the result for the single weighted 
functions is independent of whether the impurity weights are applied to the interior or to the 
exterior of the single particle resolvents, such that (Gf(l, 2)G° b i (l, 2)) = (G°J(1, 2)Gf(l, 2)). 
Indeed it is not hard to show that other arbitrary distributions of two impurity weights also 
give the same results as long as there is one weight applied to each resolvent. The same is 
of course true for host related properties. 

Accordingly, b k from ([49]) in the previous section goes over to 

S a S fe + A(a k - a k ) 
h = a k — — (72) 

and similarly c k from equation (0) to 



16 



c k = S a S fe a fc ( 1 + A-— ^-r— ) . (73) 
We remember that the Fourier transform of the unweighted two particle Greens function is 



£e^(G a (l,2)G b (l,2)> = °* . ■ (74) 



Phenomenologically, one can scrutinize the uncorrelated limit where A — > for which in 
( [73]) the fc- dependent part of c k behaves as 

hm c k = S S 6 a fc + ... (75) 

which is the correct limiting result for a product of two averaged double impurity weighted 
single particle functions 2)G"(1, 2) or the weighted two particle function without co- 

herent corrections. 

However, taking this limit is reasonable only in special cases since generally the vertex A 
depends on the self energy S. It will be appropriate to take in a weak disorder limit (virtual 
crystal limit) since for V a / b — > the CPA predicts that T, a / b — > cV a / b and A — > c(l — c)V a V b 
which means that A approaches zero faster than E a /6 in this limit. 

One should also note at this point, that although it is apparent that the term (l + S a F a + 
S5F5) , which arises in the difference of the Fourier transforms a k and a k , originates from 
the difference in the diagonal parts of their respective real space Green's functions, it can 
not be neglected in this treatment. When A is finite, this term is multiplied with other 
k— dependent quantities and thus contributes to the off diagonal elements as well. Written 
in terms of a k , ( [73] ) can be cast into 

c k = - v Qfc - (s a S 6 + ^ b (a k )-\l + E a F a + S b F b )) (S a S b + A(l + S a F a + S a F b )) . 
L a L b - Aa k v ' 

(76) 

If again diagonal terms in real space are neglected, i.e. k— independent quantities in space, 



it can be shown that ( 76 ) can be cast into the compact form 



Ck = v v Qfc . (S a S b + A(l + S a F Q + S 6 F 6 )) 2 (77) 
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where we have introduced (k = C& + c G , whereby c D is independent of According to the 
definition of this finally relates the weighted off-diagonal real space two particle Green's 
function to the unweighted one as 

<G?(1,2)G*(1,2)> = (G B (1, 2)G b (l, 2)) + A(1 + g'f° + SfcFfc))2 . (78) 

{v a v b y 

Equally, the relation for the single weighted function (G % a (l, 2)G l b (l, 2)) can be obtained 
almost immediately if ( |56l) is modified for this choice of and C^, which goes over to 

d k = a k + Aa k - " fc = - - Qfc [E a E b + A(l + E a F a + Z b F b )} (79) 

thus relating the real space functions in this case as 

(Gl(l, 2)G*(1, 2)) = (G a (l, 2)G b (l, 2)) SflS * + Ml+^g^) (80) 

where no diagonal contributions were omitted. 

The fact that in the single weighted case the same weighting factor occurs once as 
opposed to twice for the double weighted one is structurally equivalent to the results for 
the single particle theory. Although retrospectively that might not be surprising, it is also 
quite interesting in view of the fact that for the most general case of the previous section 
the weighting factors in (|50D and ( 57|) almost look alike were it not for the difference in the 



pre-factors, i.e. a k in the former and afc/E a Ef, in the latter case. 

In principle all other functions can be derived now using the set of relations (|63|) - (|65D . 
However, one can save a considerable amount of algebra by recalling the following relation 
which holds in the transition from impurity to host related properties in the CPA 

-E^(y-E) c _>(i_c). (81) 

Starting from the definition of the vertex A in ( |57|) and the CPA condition (T n ) = one 
finds that 

n -i 



A 



1-c 



(K-E a )(H-E 6 ) E a E b 



(82) 
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which in a single band case (a = b) readily simplifies to the form first introduced by Leathi 

Aa=b ~ 5F~ 1-(V-2E)F- 

With the relation (|6]) between F and £ the weighting factor for the single weighted impurity 
function from fl80|) reduces to 

S a S b + A(l + S a F a + S b F b ) = (1 - c)A 

V a V b (V a -E a )(V b -E b )- [ ' 

Correspondingly, the factor for the double weighted function is the square of this quantity. 
One can now employ (^3]) and (|64l) to find the weights for the single host- and the impurity- 
host weighted functions, respectively. We find 

rA 2 

(GiGf) = (G a G b )— (85) 

^a^b 



and 



= ^ E.E^-E^-E.) < 86 » 



which by means of ( [55]) gives 

cA ^ 2 



(Gf ) = (G a G b ) ( — ) . (87) 

^a^b' 



Note that this could also have been expected from the transformation property (plf) . The 
impurity and host weights as represented in ([84] ) and (p5|) are immediately seen to be the 
two contributions summing to A -1 in (p2|). Thus, even though there are many different 



possible representations of the two band vertex, the representation in fl8~2|) shows that the 
transition from impurity to host properties leaves A invariant. 

A notable feature about the weights calculated in this section is that they are independent 
of any wavevectors and only multiply the unweighted function as scalar energy dependent 
factors. This is a direct consequence of the single site approximation. 
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B. Theory of Absorption 



In this section we evaluate a form of the two particle functions needed in the calculation 
of the linear response absorption in a disordered solid. We take the operators C^ 1 ' and in 
the characteristic form of dipole operators similar to the one in (|B2"|). Furthermore, we assume 
that the dipole matrix elements be essentially constant, such that C^ 1 ''^ = ^^/^ Y^ m \ 
m){m |. This choice corresponds to applications for the description of processes involving 
transitions between bands of different angular momentum symmetry such as required by 
the selection rule for optical processes cit £it zero total momentumE'0. At the end of this 
sub-section we will also give for completeness a short account of the calculation of linear 
response conductivities in disordered solids. 

For the above choice of 7^)/( 2 ) the calculation of the term K, in (^) is greatly simplified 
compared to before, but it will be necessary now to also consider total diagonal terms, since 
the sums over all states in C^ 1 ^^ couple all sites and hence all contributions coming from 
term Ai in (|2~3| ) have to be included. 

As a consequence of the introduction of the dipole operator the main change arising in 
the result for K, is that the site diagonal elements (n | JC \ n) and (n \ G a C^Gb \ n) as well 
as (n | G a g~ 1 C^g b ' 1 Gb \ n) are now actually independent of nEl. Since 

C(i) 

now couples the 

functions to its left and right like a matrix product, the B-S equations (|38| ) and (^) have a 
very simple solution in terms of their Fourier transforms. Introducing a,b,c = cik,bk, Ck \k=o 
as the zero momentum elements of the respective transforms from last section we effectively 
get c = (n | K, | n) , a = S a Sb(n | G a Gb \ n), where the omission of C^ l > indicates that 
the two single particle resolvents are now simply multiplied as matrices, and equation (f49|) 
reduces to 

E a S 6 + A(a - a) 



S a S b - Aa 



(88) 



and equation (loTf ) to 



a 



t = w l 1 + feJ (89) 
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which can be recast into 

c = - ° - [S a S fe + A(l + S a F a + E b F b )} 2 
+ (1 + S a F Q + [£ £ 6 + A(l + S a F a + (90) 

where the term independent of a which was discarded in the previous section has to be kept 
in this case since the contributions of the diagonal elements become important. 

In calculating Ai one has to be aware that U is a matrix which just has a random 
occupation of its diagonal. The sum over all sites in the operators and will hence 
just pick out the sum of all total diagonal parts (n | Q | n) (n | U | n) where Q is a generalized 
product of several Green's functions in the same band (the case where U and Q are swapped 
is analogous). From (21) we get immediately 



M = cV a V b + cg- l V b - cg^Gag-W, - g-^T^g- 1 ^) + (a <- b) (91) 



where the only term giving slight complications is (T a G a g a 1 Ub). However by means of ( |19|) 
we obtain 

{T a G a g- l U b ) = Y, a (T a G a U b ) + (T a U b ). (92) 

In this expression, the second term presents more complications, since for the first one we 
remember from ([|) that 

TG = G~ X G - 1 (93) 

and Ub = U a Vh/V a such that 

(T a G a U b ) = (G- x G a U b ) -cV b =^±- cV b . (94) 
Subsequently, the second term can be decoupled by means of (|27|) - (|29|) 

(T a U b ) = 52(Q a n U b ) = + E Q a m G a )T:U b ), (95) 

n n m^n 

Applying (|30|) and using the the fact that U is diagonal yields 
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(Q a m G a TZU b ) = (Q a JG a (T:U b ) = 



(96) 



since (m ^ n) and (Q m ) — and hence we find 

(T a U b ) = £<T n ^> = - C{ya ~ ^ . (97) 

n 1 — {V a — l^aj-Ta 

Collecting all terms for M. and some more algebraic manipulation finally yields 

M = cV a V b - (1 + E a F )E a V 6 - (1 + E 6 F 6 )E 6 V a . (98) 



The sum of all diagonal parts, i.e. M. and the ones from the second term in (|90"D , can be 
shown to assume the very compact form 

AF a F b 



M + (1 + S a F a + E b F b ) [E a E b + A(l + E a F a + E b F b )] 



(99) 



The final result for (G^G l b l ) thus amounts to 

(1 - c)A 



AF a F fe 



(100) 



(K-E a )(y 6 -s 6 ) 

Here the first term has been rewritten in the same way as already derived in the last sub- 
section for the finite range susceptibility. From there it is also seen that the single weighted 
function will have the same weight as calculated in (|84l) for the corresponding function in 
the exchange coupling case. We find 

(l-c)A 



{G\G\) = (G a G b ) 



(K-E a )(H-E 6 ) 



;ioi) 



and by means of ( |63D - (|65|) 



(G a G b ) — (G a G b ) 



cA 
E a Eb 



(102) 



(GfGl h ) = {G a G b ) 



c(l - c)A 2 



AF n F b 



(K-E a )(H-E 6 )E a E 6 VaVb 



(103) 



(Gf Gf ) = (G a G b ) 



cA 



E„E/- 



2 AFaFl 
V a V b ' 



(104) 
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Comparing the results of the last two sub-sections, it becomes clear that the several 
weights obtained are essentially universal. The main difference in the absorption case as 
compared to the susceptibility one comes from the diagonal terms which have to be kept in 
the double weighted functions. The single weighted analogues are void of this difficulty and 
the weighting factors are identical for both cases. 

So far we had omitted to consider a form of the two particle functions which is needed 
for conductivity calculations. However, VelickyB showed for the unweighted functions that 
the vertex T vanishes in the corresponding expression for the conductivity, due to the anti- 
symmetry of the dipole matrix elements in /c-space if they are taken between Bloch states 
of a non interacting Fermi system in a crystal with inversion symmetry. The same is also 
true for the weighted case and effectively the weighted functions which would have to be 
used for such calculations would just consist of products of the corresponding single particle 
quantities. 

This however turns out to be a general deficiency of the CPA in the single site approxima- 
tion since due to the multiple scattering exclusions, only ladders of nested diagrams are used 
in calculating the total contribution of the coherent scattering. The CPA therefore neglects 
higher order two particle correlations which are in fact non-zero and contribute markedly 
to the conductivity. Langer and Nea@ have shown that the so called "maximally crossed" 
diagrams, i.e. diagrams which have a maximal crossing of coherent particle-particle scatter- 
ing lines actually contribute the leading part - in the order of the expansion considered - 
to the full two particle disorder vertex for the conductivity in an otherwise non-interacting 
system. For the case of interacting Fermi systems, however, the presence of the interactions 
is sufficient to destroy the aforementioned antisymmetry and thus also the terms already 
included in the vertex of a single site two particle CPA as discussed here will give a finite 
contribution to the conductivity in real systems. 



23 



C. Split band limit 



As already indicated, in contrast to many other theories of disorder, the CPA interpolates 
correctly to the limits of strong disorder and high concentrations. In this situation each band 
splits into two components of strength c and 1 — c, respectively, which represent largely 
separate A— and -B-type excitations. An intuitive consideration of the underlying physics in 
this limit suggests, that the correct description of an absorption process should predict that 
the overlap integral for transitions between sites pertaining to different alloy components 
will gradually decrease and thus in reverse, that transitions between sites of equal type 
will be more and more favored. In the following, we prove that the CPA of weighted 
two particle functions predicts this behavior correctly which makes it useful for a better 
quantitative understanding of absorptive and dispersive processes in strongly disordered 
alloys. The corresponding single particle theory must fail in this respect, since it will weight 
the occurring transitions only with the products of concentrations of sites involved in these 
transitions. 

To illustrate this we assume that our material components A and B have corresponding 
single site energies e^, £^ an d e a > e b f° r the two bands respectively and the carriers have 
become totally localized, i.e. their effective mass goes to infinity, or vice versa the bandwidths 
involved go to zero. The potentials V a and V b are then defined as 

V a = et- e B a , V b = ef- e? (105) 

and the single particle site diagonal Green's functions go over to 



c 1 — c 1 

Fx(z) = —j + —b = -g (106) 

Z — € x Z — £ x Z — £ x — l^x 



where A labels the corresponding band. Thus, if a two particle theory is constructed from 
single particle properties only, and coherent terms in the two particle scattering are neglected, 
this leads to peaks in the absorption spectrum as shown in table |. The energies at which 
the peaks are centered are shown in row 2 and their relative weight for the uncorrelated 
average is shown in the row 3. 
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In the following we show for the limit of strong disorder how, upon introducing the ver- 
tex corrections in conjunction with the appropriate weighting factors for the two respective 
components, the expected transitions are filtered out correctly with their appropriate tran- 
sition strengths and the spurious crossed terms are suppressed as shown in row 4 of table |. 
Defining 

z-ef = x (107) 

F and £ can be rewritten as 

, 1 x-(l-c)V , s 

= v = / T/ 108 

x — S x{x — V) 

CvT, 

EW = T^W^j (109) 

and the vertex A given by ( |3^) can be written 

(1 - c)cV a V b x a x b {x a - V a )(x b - V b ) 

~ [x a - K(l - c)] [x b - H(l - c)] [cx a x b + (1 - c)(x a - V a )(x b - V b )] ' 1 J 

The impurity weight £ = (1 — c)A ct) /(V c — S C )(K — S„) can be represented as 

£= (111) 

cx a x 6 + (1 - c)(x a - V^)(X6 - Vfe) 

and equally the host weight rj = cA afe /S a S fe as 

(l-c)(a; a -K)(^-H) (n2) 
cx a x 6 + (1 - c)(a; a - K)(a; 6 - H) ' 

It is now evident, that the impurity weighted quantities are proportional to c and the host 
weighted ones to 1 — c and not the other way around as their the appearance in terms of their 
weighting factors might superficially have suggested. The correction factor coming from the 
diagonal terms in the double weighted functions can be recast into 

73 *£5 = , (1 ~ c)c w ? . (113) 

V a V b CX a X b + (1 - C) {X a - V a ) {X b ~ V b ) 

Furthermore, the unweighted but vertex corrected two particle propagator K assumes the 
form 
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K = CXgX b + (1 ~ C)(gg ~ Vg)(X b ~ Vb) / 114 S 

X a Xb{x a ~ V a )(x b ~ V b ) 

After further algebra one can show, that in this limit the weighted functions can be expressed 
as 

< G « G «) = K? + , = «*?> = = (z> _ yj) C (zt _ yt) (H5) 
(G hh G hh ) = Kt 1 2 + 1 = (G h G h ) = K V = (116) 

which is exactly what is expected to happen physically in this limit. The crossed terms 
in the transition are canceled out - hence the crossed function in ( P-17| ) goes to zero - 



and the double weighted functions become identical to the single weighted ones, since now 
effectively only the totally site diagonal element Kp = (G utJ, (l,l)G UfJ- (l,l)) still contributes 
to the transitions, which implies that only two possibilities for weighting the two particle 
functions remain, namely as (G a G l b ) and (G^G b ). As could be expected from a theory which 
properly describes the strong disorder limit the transition strengths now distribute with the 
concentrations c and 1 — c between the Ab — > A a and Bb — > B a transitions, respectively, 
such as shown in row 4 of table |. This feature may in reverse be used to derive the total 
diagonal element for all ranges of disorder. By requiring K l £ = we find 

Kd = x = ggfg i * + P'.-sjw-*)V 

(V VM \ c (1-c) ) 1 ' 

The weighted versions of this element are obtained by just multiplying the corresponding 
single weights from (|101| ) and (|102| ) on to it. Moreover, Kr, is equivalent to the = 



component of the two particle function calculated for the finite range susceptibility in the 
last sub-section. In terms of the notation introduced there it reads 

which would have been harder to evaluate starting from that representation. The total site 
diagonal element thus decouples into the corresponding site diagonal single particle functions 
with an appropriate correction term accounting for the coherent processes. 
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V. NUMERICAL RESULTS 



In this subsection, in order to exemplify the general results, we discuss numerical results 
obtained for the optical absorption in a non-degenerate binary semiconductor alloy for a 
given model density of states. We are thus able to show how a CPA type of polarization, 
including vertex corrections, decomposes into contributions originating from single alloy 
components as the strength of the disorder is increased thus eventually causing the joint 
density of states to split into several components (up to three different ones for the double 
weighted case). 

In all our calculations we have used a semi elliptic density of states for a pair of 3-D 
conduction and valence bands as introduced in Ref. |6] for the single particle CPA, i.e. 



where A labels either the conduction A = a or valence A = b, and w\ is the half-width of the 
band considered. This has the advantage that the self consistent CPA equation for the self 
energy T,(E) is a third degree polynomial which can be solved analytically. 

To understand the effects that arise from genuine two particle behavior as compared 
to those expected from the single particle CPA, we recapitulate some of the features of 
the single particle theory first, mainly building on the treatment presented in Ref. || It is 
established there, that depending on the concentration and disorder strength relative to the 
bandwidth, an impurity band is eventually split off while in this split regime under some 
circumstances the CPA self energy exhibits a pole. Fig. [l| shows a reproduction of the 
"phase" diagram first presented there, indicating how the several regions are separated. It 
can be seen that for a disorder strength | V\/w\ |> 1 the bands always split into A— and 
S-components, whereas the splitting occurs earlier as the concentration c is reduced, going 
down to | V\/w\ |> 0.5 as c — > 0. 




(120) 



Px(E) 







E \<w x 



27 



We have calculated the linear polarizability of the medium by employing a Kubo formula 
as introduced in fl59|). Furthermore, we continue assume that the optical matrix elements 
are essentially constant and that such elements are the same for both alloy components and 
we hence normalize them to unity. The optical absorption is the negative imaginary part 
of the retarded polarizability of the disordered medium — ImII(a;), which can be formally 
written as 

/d k 
j^iGaik; ie)G b (-k; iu - ie)) (121) 
(2tt)^ 

whereby the /c-integration is understood to be carried out after the configurational average 
has been performed, since before that both resolvents would depend non-trivially on two 
momentum variables. 



We consider our system at zero temperature and follow partly the method used in Ref. |18 



for our calculations. At T = 0, the polarization can be obtained as the energy convolution 
around the conduction band branch cut of the fc-summed vertex corrected two particle 

/d^ k 
- — —(G a (k]Zi)G b (—k]Z2)) ) such that 

Ur(u) = <f K(z, uj + i5- z)dz (122) 



where we have taken over the following definitions from Ref. 18: 



ts( \ R{zi,z 2 ) , s 

1-A( SI , (123) 

where A(^ 1; z 2 ) is the usual CPA vertex from ( [82]) and 

/d h" — — 

—G a (k, Zl )G b (-k,z 2 ) (124) 

is the average-decoupled two particle function. Assuming that the conduction and valence 
band dispersion relations exhibit a similar shape such that they scale proportionally 

e_M = T £ J^1 (125 ) 

W a W b 



( |124j) can be shown to simplify to 
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, w a F a {z 1 )±w b F b {z 2 ) 

K{Zi,Z 2 - ; —, — rr— ; — ; — rr {L*0) 

W a [Z 2 - ^b{Z2)\ ± w b [zi - S a (Zi)J 

where F\(z) are the site diagonal single particle functions first introduced in connection with 
(p|). As usual we assume that the effective mass of an electron in the conduction band is 
positive and that of a hole in the valence band is negative. Accordingly we have chosen the 
upper choice of signs in ( [L25| ) and ( |126| ) for our calculations. 



To be able to analyze the obtained results with regard to the effect of the the inclusion of 
vertex corrections we first consider qualitatively the features that would be expected from the 
transition process represented by the energy convolution in (|122[) in an intermediate regime 
of disorder, if the configurational average in the two particle function is decoupled and 
effectively only single particle properties are employed. This would correspond to replacing 
K from (|123|) by R from (|124|) in ( |122j ). We assume for now that the concentration be about 



0.5 and the bands have just split by a notable amount. With the semi-elliptic bands used, 
the transition process can be represented as shown in Fig. [2| 

The disorder strengths give approximately the distance between the centers of the single 
bands. The convolution of two separated finite bands, occurring in R, would yield a set of 
finite bands in the joint density of states (DOS) whose width is the sum of the widths of 
the contributing components. Two cases are considered where the band offsets of the A— 
and 5-components of the alloy are in equal or opposite directions corresponding to parallel 
or anti-parallel disorder. 

In the case of parallel disorder, this would amount to the A b — > A a and B b — > B a 
-transitions lying in the center of the joint DOS, framed by the contributions from the 
B b — > A a and A b — > B a -transitions upper and lower end respectively as shown in Fig. || 
At c = 0.5 these regions would have relative distribution of weighted states of 1 : 2 : 1 from 
lower end : center : higher end. In the case of anti-parallel disorder the picture should be 
similar with the only difference that the spectrum is turned inside-out with the B b — > B a 
and A b — > A b components on the top and the bottom end of the joint DOS and the mixed 
transitions in the center, again with a distribution of 1 : 2 : 1. 
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The calculations in the previous section for the split band case strongly suggest that the 
vertex corrections will increasingly suppress the cross transitions as the disorder strength 
is increased, which is verified in our numerical results. Indeed, our results show that this 
suppression is already displayed quite strongly in an intermediate disorder range, i.e. in a 
regime where the single bands just begin to split. 

The appended plots in Figs. £| and || show cumulative absorption spectra calculated 
from (|122|) as well as their single and double weighted components for parallel and anti- 
parallel disorder of various strengths, covering both the joint and the split band regime. 
The conduction band half-width is normalized to unity and the valence band half-width 
is taken to be 0.8. The concentration of impurities is fixed to 0.35 in order to study the 
high-concentration behavior rather than dilute impurity effects. 

In the transitional region when the disorder strengths start to exceed the single particle 
half-bandwidths | V\/w\ |> 1 and the conduction and valence bands start to split we observe 
the following behavior: In the case of parallel disorder shown in Fig. |] the spectrum starts 
to exhibit a discontinuity in its derivatives at the flanks accounting for a pair of mixed 
components splitting off sideways from the main contribution. At the transition point the 
contributions of the flanks relative to the central bulk part are approximately distributed in 
a ratio of 1 : 14 : 1. In the case of anti-parallel disorder shown in Fig. [| the suppression of 
the mixed transitions is even stronger so that as soon as the single bands split the crossed 
transitions cancel out entirely within the numerical accuracy and the joint density of states 
starts to exhibit a gap between two separate contributions which mainly consist of A— and 
^-transitions, respectively. 

In a regime of strong disorder, further beyond the splitting of the single bands as shown 
in the last plots of Figs. f| and |5], one finds that over large regions the single and uniformly 
double weighted contributions coincide almost exactly, implying that the spectrum is built 
almost entirely from the total diagonal element of the two particle function, which we had 
calculated at the end of section IV. The total diagonal element can be obtained as an 
independent /c-sum over the two single particle resolvents involved. This situation represents 
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a breakdown of the fc-selection rule which holds in pure media. 

If the splitting of II into single weighted components II = II + U B is considered and 
the result is compared in appearance with the splitting into components of the site diagonal 
single particle functioni 



one finds that the single and the two particle behavior appear to be strikingly similar in the 
case of anti-parallel disorder, as can be seen if the plots for the imaginary part of the single 
particle function shown in Figs. |6] (a) and (b) are compared to the ones for the corresponding 



In a system without disorder this similarity is evident if there are no further local inter- 
actions in the problem, since the non-interacting two particle motion decouples into a center 
of mass and a relative co-ordinate and while the center of mass motion can be set to zero, 
the relative one can be mapped onto a single particle co-ordinate. Upon the addition of the 
disorder, this decoupling fails to work and it can only be regained by using an appropriate 
configurational averaging procedure. However, if in the presence of disorder an average is 
only performed on a single particle level, thus omitting average induced two particle cor- 
relations the reduction obviously fails to work, as is seen through comparing the plots of 
Fig. [7] showing a spectrum calculated without the vertex corrections with the ones for the 
corresponding parameter values of Fig. |5], which properly include these corrections. The 
results show, that in regimes of intermediate and strong disorder the influence of the vertex 
corrections is very substantial. 

In comparison to the single weighted (two fold) splitting the double weighted (three 
fold) splitting exhibits a rather curious behavior. Even though the components — Imll 
and — lmU BB lie underneath their single weighted complements -Imll' 4 and —lmU B in 
some parts of the spectra, which one might expect to happen globally on first thought, 
they either coincide with them or even exceed them in other parts - sometimes to such an 
extent that they reach beyond the cumulative function. However, it has to be noted, that 




(127) 



parameters in the two particle case of Fig. |5] (third and fourth plot from the front). 
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these components, like the unweighted function, are always uniformly positive in sign, and 
therefore exhibit the correct analytic behavior that a function defined on this footing has to 
satisfy. It is required that these components be positive definite, because the net absorption 
in the medium must always be positive unless the system is excited out of equilibrium, which 
we do not consider here. In the preceding discussion we have assumed for convenience that 
the optical matrix elements between states of different components of the alloy are equal. If 
different optical cross sections (matrix elements) fi A ^ B are distinguished for A— and 5-atoms 



the weighted contributions to yield the integrand of (|122|) would sum as 

K = (fi A ) 2 K AA + 2fi A fi B K AB + (fi B ) 2 K BB . (128) 

This shows that it would be possible to observe one of the functions U AA or U BB predom- 
inantly if either fx A or fi B happens to be much larger than the other. The mixed function 
U AB , however, will never be a separately observable quantity in a general case, no matter 
how the cross sections fi A and fi B scale relatively and therefore not so rigid criteria for its 
analyticity apply as for the uniformly double weighted functions. 



VI. DISCUSSION AND CONCLUSION 

In the previous sections we have obtained expressions for a wide class of weighted two 
particle Green's functions. The large choice for possible weightings is substantially reduced 
as restrictions are made to functions which would be useful in linear response theory. In both 
cases which are discussed for this kind of application, only five different weighted functions 
remain of which only two are genuinely independent. 

The structure of the weighting process is closely related to that derived for the single 
particle theory with the main difference that now the weights also depend significantly on 
the CPA vertex corrections. The calculation for the split band limit, the domain in which 
the CPA is superior to most of the other theories of disordered systems, gives a direct insight 
into how the properly weighted CPA extracts the correct limiting behavior from different 
possible physical processes. 
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Some care is needed in interpreting the precise physical meaning of the weights, since 
they are obtained for the averaged functions, which describe the disordered medium as 
effectively homogeneous. The concept of the propagation of a particle between sites of 
different components is therefore lost in the effective medium as a consequence of averaging 
and the initial exclusion of specific propagation paths in the unaveraged function leads to 
the effective weights. These weights simply account for the average partition in probability 
for the simultaneous propagation of the particles between partly or completely specified site 
types at a given pair of energies z\ and z 2 . 

Our numerical results show the general importance of the inclusion of vertex corrections 
into a properly self consistent two particle formalism. We managed to visualize that, as a 
consequence of the inclusion of these average induced two particle correlations, the center of 
mass and relative motion of the two particle system effectively decouple to a large extent. 

We believe that the general method developed here will find applications in various 
situations where two particle motion is studied in a disordered medium. The effect of alloying 
on the electronic susceptibility and hence for example on the Rudermann-Kittel interaction 
has already been mentioned. In particular we believe it can be extended for use in systems 
where the two particles have a direct interaction, such as the Coulomb interaction between 
carriers occurring in excitons in alloyed systems. In such a case the static correlations 
between particles, created by the disorder and accounted for by the vertex corrections, and 
the dynamic correlations introduced through the carrier-carrier interaction, create additional 
static-dynamic correlations. Moreover it may be possible that the underlying disorder of the 
system gives rise to an induced disorder to the carrier-carrier interaction itself. Both of 
these effects can be treated within the method developed here and will be addressed in a 
forthcoming publication. 
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TABLES 
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TABLE I. Energies of poles at which interband transitions may occur versus center of gravity 
weights in an asymptotically strong disorder limit before and after the inclusion of vertex correc- 
tions. 
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FIGURES 

FIG. 1. Regimes for the CPA of the single particle density of states (120) for w\ = 1 in a 
binary substitutional alloy depending on the impurity concentration c and the disorder strength 
V. In the lowest region the band is unsplit. Above the first broken line the bands split into two 
components and above the second broken line the CPA self energy S exhibits a pole between these 
components. 

FIG. 2. Diagrammatic representation of the energy convolution in (122) effectively using the 
term R(z±, z-i) from (126) only. If R is taken at Z\ = Z2 and c is about c = 0.5 the conduction and 
valence bands are split and the imaginary part of R decomposes into four components centered 
about e„ , and as shown qualitatively. The components are separated by the band gap 

Eg in the middle and the single particle splittings in the conduction and valence parts at the upper 
and lower end, respectively. The energetic difference between the gravity centers are given by 
V a = £q — and Vb = e£ — ef , as already introduced, as well as A = — ■ Bracketed 
expressions denote anti-parallel disorder. 

FIG. 3. Qualitative result of the convolution (122) as it can be expected if performed using the 
term R(zi,Z2) from (126) only. At about c = 0.5 the joint density of states would be expected to 
distribute over regions in terms of their compositional origin as shown. The inclusion of the vertex 
is expected to primarily suppress the mixed ^45-transitions and favor the AA— and i3S-transitions 
as the disorder strength is increased. For parallel disorder this leads to a suppression of the spectra 
at the flanks while for anti-parallel disorder (in brackets) it causes a suppression of the central 
part. Both these features are well represented in the calculated spectra of Figs. 4 and 5. 



37 



FIG. 4. Negative imaginary part of the polarizability of the disordered medium (joint DOS) for 
parallel disorder and its possible weightings into A— and B— or AA-, AB— and Ui?-components. 
The concentration of ^4-atoms is fixed to c = 0.35 and the band half-widths were taken to be 
w a = 1.0 and Wb = 0.8 for the conduction and valence bands, respectively. The disorder was varied 
through a range of parameters, annotated on the bottom right of the plots, with particular focus 
on the splitting point | V\/w\ |= 1 of the single particle bands. At the splitting point a pair of 
flank-components starts separating of sideways and gets completely isolated as V\/w\ — > 3 in the 
last plot. These flanks are constituted to about 50% by the AI?-component, whereas the central 
region almost entirely consists of AA— and l?i?-components only. 

FIG. 5. Negative imaginary part of the polarizability for anti-parallel disorder (reversed in the 
valence band) with otherwise equal model parameters to the parallel case shown in Fig. 4. The 
joint part starts to exhibit a gap as soon as the conduction and valence bands do. In the split 
regime the spectrum almost entirely builds up from single A— and ^-components which in the 
interior of the two bands again coincide largely with the AA— and SB-components indicating a 
strong breaking of the /c-selection rule already in an intermediate disorder regime. In the last plot 
V a = 3.0, Vt, = —2.4 it can be seen that the double weighted functions are non-zero in the gap 
region. This means that in a case where the optical cross sections \x A l B of the of A— and S-atoms in 
(128) are non-equal, the gap in the observed absorption spectrum would be less complete containing 
either AA— or SB-states primarily. 

FIG. 6. Imaginary parts of the site diagonal single particle function F taken for the same 
concentration c = 0.35, a band half-width of w = 1.0 and disorder strengths of (a): V = 0.8 and 
(b): V = 1.2. By comparing the splitting behavior of this function with the one of the corresponding 
two particle function for anti-parallel disorder, shown in the third and fourth plots in Fig. 5, one 
finds that the behavior, particularly the one of the single weighted A— and B-components, is 
strikingly similar. 
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FIG. 7. Negative imaginary part of the polar izability without vertex corrections for anti-parallel 
disorder. The parameters taken were (a): c = 0.35, V a = 0.8, = —0.64 and (b): c = 0.35, 
V a = 1.2, Vb = 1.0. The characteristic splitting of the joint DOS, which is displayed in the properly 
corrected function of Fig. 5 as the disorder strengths exceed the value | V\/w\ |= 1, does not 
occur and the shape of the joint DOS rather resembles the qualitative one of Fig. 3. The curves 
for the single weighted components in both (a) and (b), which were obtained from correspondingly 
weighted uncoupled single particle functions, do not satisfy the sum rule II = IT 4 + IT 3 . 
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